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Abstract: Investigation of the degradation kinetics of perfluorooctanoic acid (PFOA) 
has been carried out to calculate rate constants of the main elementary reactions using 
the multichannel Rice-Ramsperger-Kassel-Marcus theory and canonical variational 
transition state theory with small-curvature tunneling correction over a temperature range 
of 200-500 K. The Arrhenius equations of rate constants of elementary reactions are fitted. 
The decarboxylation is role step in the degradation mechanism of PFOA. For the 
perfluorinated carboxylic acids from perfluorooctanoic acid to trifluoroacetic acid, 
the quantitative structure-activity relationship of the decarboxylation was analyzed with 
the genetic function approximation method and the structure-activity model was constructed. 
The main parameters governing rate constants of the decarboxylation reaction from the 
eight-carbon chain to the two-carbon chain were obtained. As the structure-activity model 
shows, the bond length and energy of C1-C2 (i?ci-C2 and Eci-ci) are positively correlated to 
rate constants, while the volume (V), the energy difference between £homo and -Elumo (Aii), 
and the net atomic charges on atom C 2 (Qcz) are negatively correlated. 

Keywords: electrochemical degradation; PFOA; rate constants; quantitative 
structure-activity relationship 
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1. Introduction 

Perfluorinated carboxylic acids (PFCAs) have been widely used in industry as surfactants, surface 
treatment agents, metal coating, fire retardants and carpet cleaners for many years [1-4]. As a typical 
perfluorinated acid, the production of perfluorooctanoic acid (PFOA) worldwide exceeded 1000 metric 
tons in 2004 [5]. PFOA is a new persistent organic pollutant and difficult to decompose in the 
environment [6]. Due to its long-range oceanic and atmospheric transportation ability, PFOA has been 
detected in various environmental matrices such as water, dust, sediments and living tissues [7-9]. Research 
shows that levels of PFOA in wildlife range from 0.05 ng/mL in the blood of cod collected from 
European waters to 8.14 ng/mL in the plasma of loggerhead sea turtles from North America [10,11]. 
It can be also accumulated in creatures through water, food and atmosphere, causing the decrease in 
fertility rate, birth weight and other immune system diseases [12-18]. Although the environmental 
protection agency of the United States (US-EPA) and the European environment agency (EEA) have 
adopted an industrial program in order to reduce the global emissions of PFOA [19,20], the remained PFOA 
in the environment still has potential risk on humans. So, it is important and urgent to find an effective 
degradation method. 

Recently, electrochemical treatment of PFOA has been considered one of the most promising 
methods due to its strong oxidation and environmental compatibility [6]. What is more, Hoffman et al. 
found that PFOA-contaminated drinking water is a significant contributor to PFOA levels in serum [21]. 
Researchers have found that PFOA could be degraded over boron-doped diamond (BDD) film electrode 
and Ti/Sn0 2 -Sb-Bi electrode [22-24]. The reaction mechanism is proposed as the direct electrochemical 
oxidation cleaves the C-C bond between the C 7 Fi 5 and COOH in PFOA to generate a C 7 Fi 5 radical and 
CO2 firstly [24]. Zhuo et al. studied the electrochemical oxidation of PFOA using Ti/Sn02-Sb-Bi anode, 
and they found that after 2 h electrolysis, over 99% of PFOA was degraded with a first-order kinetic 
constant of 1.93 h~ . Then, the degradation mechanism was revealed according to the intermediate 
products detected [23]. Lin et al. studied the electrochemical degradation efficiency of PFOA using 
different anode, current density, pH value, plate distance and concentration. They also put forward 
a reaction mechanism in accordance with experiment results [6]. Although much work has been done to 
study the degradation effect of electrochemistry, the degradation mechanism has been partly revealed, 
and the rate constants of elementary reactions have not been reported so far. The reaction rate is helpful 
to find an optimal reaction way, then a kinetics study is needed. Furthermore, the quantitative 
structure-activity relationship (QSAR) analysis is helpful to understand the reaction mechanism. 

In this study, the geometrical parameters are optimized at the MPWB1K/6-31 + G(d, p) level. On the 
basis of the quantum chemical information, the rate constants are calculated using the multichannel 
Rice-Ramsperger-Kassel-Marcus (RRKM) theory and canonical variational transition state theory (CVT) 
with small-curvature tunneling (SCT) correction over a wide temperature range of 200-500 K. 
Then, the kinetics study has been performed to calculate the rate constants of elementary reactions, 
and the quantitative structure-activity relationship from perfluorooctanoic acid to trifluoroacetic acid 
with rate constants is analyzed in order to explore the main factors affecting the rate constant of 
decarboxylation reaction. 
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2. Results and Discussion 

The main possible reaction paths of electrochemical degradation of PFOA are drawn in Figure 1 . 
Figure 2 shows two reaction circles in the electrochemical degradation pathways. The chemical 
structures of transition states in electrochemical degradation reactions are shown in Figure S 1 . 

2.1. Reaction Mechanism 

The detailed electrochemical mineralization mechanism of PFOA has been analyzed in previous 
study [25]. The first step of electrochemical degradation is the electron transfer process from carboxylic 
acid to anode in which the intermediate IM1 is generated [24]. Then, perfluoroheptyl radical IM2 and 
C0 2 are produced from IM1 via the transition state TS1, which is a decarboxylation reaction. In the 
subsequent reactions, IM2 may continue to react with OH, O2 and H 2 0 in the electrolytic cell. Channel A 
shows the reactions initiated by OH. The addition reaction of IM2 with OH generates adduct IM3, which 
is a barrierless process. Perfluoroheptanoyl fluoride IM4 is formed through a HF desorption process in 
the four-membered ring of IM3 via the transition state TS2. Obviously, IM4 can be hydrolyzed to 
perfluoroheptanoic acid (PFHA) PI and hydrofluoric acid HF. What is more, IM4 also can continue to 
react with OH to generate adduct IM5 via the transition state TS3. Then, the fluorine atom will be 
removed to form PI via the transition state TS4. A HF desorption process in the four-membered ring of 
IM5 can also occur to generate PFHA radical P2 via the transition state TS5. Obviously, the HF 
desorption process of IM3 is challenged due to its high potential barrier of 53.43 kcal-mol . An easier 
way lies in the hydrogen atom of the IM3 being abstracted by OH to form intermediate IM6 via the 
transition state TS6, and the potential barrier is 6.62 kcal-mol -1 . Then, IM6 can remove a fluorine atom 
to produce the IM4 via the transition state TS7, or be decomposed into perfluorohexyl radical IM7 and 
difluorophosgene via the transition state TS8. The above reaction barriers are 33.02 and 8.32 kcal-mol -1 , 
respectively. Difluorophosgene can be hydrolyzed to the final products HF and CO2. Obviously, 
the latter path occurs much more easily thermodynamically. 

Channel B depicts the reactions initiated by oxygen. IM2 reacts with O2 to generate adduct IM8, 
which is a barrierless reaction. The oxygen atom of IM8 can be abstracted by OH and 0 2 to form 
IM9 and IM6 via the transition state TS9 and TS10, respectively. IM9 can also be decomposed into 
IM6 and OH, but considering the high potential barrier of TS9, it is a difficult reaction. 

The reaction of IM2 and H 2 0 is shown in Channel C. The hydrogen atom of H 2 0 is abstracted by IM2 
and IM10 is generated via the transition state TS11. The potential barrier and endothermic energy are 
13.23 and 3.62 kcal-mol -1 , respectively. 

As shown in Figure 2, the electrochemical degradation of PFOA has three circles. The first one is the 
degradation from PFOA radical to PFHpA radical, the subsequent degradation process is from PFHpA 
to PFHxA radical, until the perfluorinated acetic acid finally decomposes into CO2 and HF. The second 
one is the degradation from PFOA to PFHpA, and the final products are same as for the first one. 
The last one is the degradation from C7F15 radical to C6F13 radical, and C6F13 radical to C5F11 radical, 
until it decomposes to CF 3 radical, which represents complete degradation. 
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Figure 1. The main possible reaction paths of electrochemical degradation of PFOA embedded with the potential barriers Ais b (kcafmol l ) and 
reaction heats AE r (kcal-moi -1 ). 
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Figure 2. The three reaction circles in the electrochemical degradation pathways. 
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2.2. Rate Constants 



The RRKM, TST and the CVT with SCT correction method are used to calculate the rate constants. 
The RRKM method is employed to calculate the rate constants of elementary reactions without barriers, 
such as elementary reactions (3), (11), and (13). Except for the above elementary reactions of PFOA, 
the rate constants of decarboxylation reactions of PFCAs radicals (-C8-C2) are also calculated for 
PFOA when it is degraded successfully after decarboxylation reaction in experiments [25,26]. The rate 
constants of elementary reactions (2)-(15) and decarboxylation reactions (-C8-C2) at the temperature 
range of 200-500 K are reported and listed in Tables 1 and 2, respectively. The reaction pressure is 
adopted at atmospheric pressure. 



Table 1. The rate constants of decarboxylation reactions (~Cs-C2) at the temperature range 
of200-500K(CVT/SCT). 



r(K) 


k&c-ic a 


7 a 
#7C-6C 


^6C-5C " 


7 a 
fsC-4C 


7 a 
K4C-3C 


7 a 
#3C-2C 


7 a 
f2C-lC 


200 


1.13 x 


10 12 


4.21 x 10 10 


1.34 x 


10 11 


2.05 x 10 11 


3.81 x 10 11 


5.17 x 10 11 


7.62 x 10 u 


220 


1.26 x 


10 12 


5.01 x 10 10 


1.54 x 


10 11 


2.50 x 10 11 


4.25 x 10 11 


5.22 x 10 11 


7.58 x 10 n 


240 


1.39 x 


10 12 


5.80 x 10 10 


1.74 x 


10 11 


2.95 x 10 11 


4.67 x 10 11 


5.28 x 10 u 


7.55 x 10 11 


260 


1.52 x 


10 12 


6.57 x 10 10 


1.92 x 


10 11 


3.40 x 10 11 


5.07 x 10 11 


5.33 x 10 u 


7.17 x 10 11 


280 


1.64 x 


10 12 


7.32 x 10 10 


2.10 x 


10 11 


3.85 x 10" 


5.45 x 10" 


5.38 x 10 11 


7.19 x 10 11 


298.15 


1.75 x 


10 12 


7.97 x 10 10 


2.25 x 


10 11 


4.25 x 10 11 


5.79 x 10 11 


5.42 x 10 11 


7.21 x io" 


320 


1.87 x 


10 12 


8.72 x 10 10 


2.43 x 


10 11 


4.73 x 10 11 


6.17 x 10 11 


5.48 x 10 11 


7.25 x 10 11 


340 


1.97 x 


10 12 


9.38 x 10 10 


2.58 x 


10 11 


5.16 x 10 11 


6.50 x 10 11 


5.52 x 10 u 


7.29 x 10 11 


360 


2.08 x 


10 12 


1.00 x 10 11 


2.73 x 


10 11 


5.58 x 10 11 


6.82 x 10 11 


5.56 x 10 11 


7.33 x 10" 


380 


2.17 x 


10 12 


1.06 x 10 11 


2.86 x 


10 11 


5.98 x 10 11 


7.13 x 10 11 


5.60 x 10" 


7.37 x 10" 


400 


2.27 x 


10 12 


1.12 x 10 11 


2.99 x 


10 11 


6.38 x 10 11 


7.42 x 10 11 


5.64 x 10 u 


7.41 x 10 11 


450 


2.48 x 


10 12 


1.25 x 10 11 


3.29 x 


10 11 


7.32 x 10 11 


3.19 x 10 11 


5.73 x 10 u 


7.52 x 10 11 


500 


2.68 x 


10 12 


1.36 x 10 11 


3.55 x 


10 11 


8.19 x 10 11 


3.38 x 10 11 


5.82 x 10 11 


7.82 x 10" 


a The unit of the rate constant is s . 
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Table 2. The rate constants of elementary reactions (2)— (15) at the temperature range of 200-500 K (CVT/SCT). 



r(K) 




A(3) b 




k(4)° 




A(5) b 




A(6) a 


k(l)" 


A(8) b 




k(9)» 


A(10) " 


£(ll) b 




k(12) b 


A(13) " 




k(U) " 


*(15) b 




200 


1.13 x 


10 12 


3.21 x 


10 




2.80 x 


10 




1.06 x 


10 




1.64 x 10~ 20 


1.54 x 10~ 32 


1.15 x 


10" 




4.21 x 10~ 24 


9.70 x 10 6 


9.92 x 


10 




1.21 x 


10~ 58 


4.75 x 


10 




2.60 x 


io- 70 


3.73 x 10 




220 


1.26 x 


10 12 


4.33 x 


10 




1.98 x 


10 




3.71 x 


10 




4.83 x 10~ 18 


1.18 x 10~ 28 


6.70 x 


10 




7.50 x 10~ 21 


2.99 x 10 7 


1.75 x 


10 




3.50 x 


io- 56 


1.91 x 


10 




1.80 x 


10 -65 


3.63 x 10 




240 


1.39 x 


10 12 


5.58 x 


10 




2.18 x 


10 




1.06 x 


10 




1.12 x 10~ 15 


2.04 x 10~ 25 


2.96 x 


10 




3.86 x 10~ 18 


7.70 x 10 7 


2.82 x 


10 




3.85 x 


10~ 54 


6.17 x 


10 




1.98 x 


IO -61 


1.67 x 10 




260 


1.52 x 


10 12 


6.93 x 


10 




5.74 x 


10 




2.59 x 


10 




1.13 x 10~ 13 


1.12 x 10~ 22 


1.05 x 


10 




7.59 x 10~ 1<s 


1.72 x 10 8 


4.28 x 


10 




2.03 x 


10~ 52 


1.69 x 


10 




5.27 x 


IO" 58 


4.33 x 10 




280 


1.64 x 


10 12 


8.35 x 


10 




4.90 x 


10 




5.62 x 


10 




5.89 x 10~ 12 


2.49 x 10~ 20 


3.14 x 


10 




7.02 x 10~ 14 


3.46 x 10 8 


6.18 x 


10 




6.00 x 


10~ 51 


4.07 x 


10 




4.59 x 


10~ 55 


7.15 x 10 




298.15 


1.75 x 


10 12 


9.69 x 


10 




1.02 x 


10 




1.04 x 


10 




1.35 x 10-'° 


1.79 x \Q- m 


7.52 x 


10 




2.52 x 10~ 12 


6.01 x 10 8 


8.28 x 


10 




8.70 x 


io- 50 


8.24 x 


10 




9.84 x 


io- 53 


6.64 x 10 




320 


1.87 x 


10 12 


1.13 x 


10 




2.84 x 


10" 




2.01 x 


10" 




3.63 x 10" 9 


1.62 x lO" 16 


1.91 x 


10 




1.09 x 10" 10 


1.08 x 10" 


1.14 x 


10 




1.45 x 


io- 48 


1.75 x 


10 


-18 


2.84 x 


io- 50 


7.02 x 10 




340 


1.97 x 


10 12 


1.29 x 


10 


-11 


5.84 x 


10 


-20 


3.44 x 


10 


-20 


5.12 x 10~ 8 


6.01 x 10~ 15 


4.05 x 


10 


-21 


2.24 x lO" 9 


1.73 x 10" 


1.48 x 


10 


■16 


1.37 x 


io- 47 


3.25 x 


10 


-18 


2.70 x 


10 -48 


4.71 x 10 


-29 


360 


2.08 x 


10 12 


1.44 x 


10 


1 1 


3.25 x 


10 


-18 


5.56 x 


10 


-20 


5.37 x 10~ 7 


1.49 x 10~ 13 


7.95 x 


10 


-21 


3.29 x 10~ 8 


2.65 x 10" 


1.87 x 


10 


■16 


1.01 x 


1Q -46 


5.66 x 


10 


-18 


1.56 x 


10 -46 


2.58 x 10 


-28 


380 


2.17 x 


10 12 


1.60 x 


10 


1 1 


1.18 x 


10 


-16 


8.60 x 


10 


-20 


4.40 x 10~ 6 


2.64 x 10~ 12 


1.46 x 


10 


-20 


3.62 x 10~ 7 


3.87 x 10" 


2.32 x 


10 


■16 


6.05 x 


jQ-46 


9.40 x 


10 


-18 


5.93 x 


IO" 45 


1.19 x 10 


27 


400 


2.27 x 


10 12 


1.75 x 


10" 


1 1 


3.01 x 


10 


-15 


1.06 x 


10 


-19 


2.93 x 10~ 5 


3.50 x 10"" 


2.55 x 


10 


-20 


3.14 x 10"" 


5.45 x 10" 


2.84 x 


10 


■16 


3.01 x 


io- 45 


1.49 x 


10 


-17 


1.57 x 


io- 43 


4.77 x 10" 


-27 


450 


2.48 x 


10 12 


2.05 x 


10 


1 1 


2.79 x 


10" 


-12 


2.44 x 


10" 


-19 


1.59 x 10~ 3 


8.21 x 10" 9 


8.39 x 


10 


-20 


2.95 x 10~ 4 


1.13 x 10 10 


4.42 x 


10 


■16 


8.93 x 


10 -44 


4.08 x 


10 


-17 


1.64 x 


10 -40 


9.11 x 10 


-26 


500 


2.68 x 


10 12 


2.30 x 


10 




6.58 x 


10 


-10 


4.89 x 


10 


-19 


3.87 x 10~ 2 


6.45 x 10~ 7 


2.23 x 


10 


19 


1.15 x 10~ 2 


2.03 x 10 10 


6.42 x 


10 


■16 


1.35 x 


io- 42 


9.42 x 


10 


-17 


4.37 x 


io- 38 


1.00 x 10 


-24 



The unit of the rate constant is s b The unit of the rate constant is cm 3 molecule 1 s 
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For the purpose of comparison, the elementary reaction (2) is taken as an example, and the TST rate 
constants, the CVT rate constants with the ZCT and the SCT correction are listed in Table SI. The TST 
rate constant at 200 K is 3.10x10 14 s -1 , which is 284 times that of the CVT rate constant at the same 

12 1 

temperature, 1.09 x 10 s . It suggests that the variational effect is significant in this reaction, and the 
higher the temperature is, the weaker the variational effect is. What is more, the tunneling effects are also 
taken into account to compute the rate constants. It is clear that the CVT constants have no significant 
difference with the CVT/ZCT and CVT/SCT ones over the temperature range of 200-500 K. It can be 
seen that the tunneling effect plays a less important role in this rate constant calculation. According to 
the previous study, the CVT/SCT rate constants are in good agreement with experimental values in a 
large temperature range [27-29]. Then, the results of the CVT/SCT method at 298.15 K are chosen for 
discussion in this paper. The TST, CVT with the ZCT or the SCT correction rate constants of the 
main elementary reactions (3)— (15) are listed in Tables S2-S11. Due to the absence of the available 
experimental values, it is difficult to make a direct comparison of the calculated CVT/SCT rate constants 
with the experimental values for all the elementary reactions. We hope that our CVT/SCT calculations 
may provide a good estimation. 

12 — 1 

The rate constant of the decarboxylation reaction (2) is 1.75 x 10 s . Obviously, the decarboxylation 
reaction occurs quite easily after the electron transfer process due to the low potential barrier. 
The rate constants of elementary reactions (3), (11) and (15) are 9.69 x 10" 12 , 8.28 x 10" 17 and 
6.64 x 10~ 31 cm 3 • molecule -1 -s" 1 , respectively, which means that Channel A is the main reaction pathway 
since the rate constant is much higher than those of the other two pathways. As Figure 1 shows, the 
subsequent reactions of the IM3 in Channel A are divided into two paths, (4) and (8), the rate constant of 

24 i — 22 3 — 1 1 

reaction (4) is 1.02 x 10 s , while the rate constant of reaction (8) is 7.52 x 10 cm -molecule -s . 
The latter is 737 times higher than that of the reaction (4), so the IM3 is easier to react with OH to form 
the IM6. Then, reaction (10) occurs readily with the rate constant of 6.01 x 10 s , which demonstrates 
that the circle of QJ^+i — > C„-iF 2 «-i is the optimal reaction pathway in the degradation process due to 
the high rate constants. The circle of C„F 2 „+iCOO — > C„-iF 2 „-iCOO is achieved by the reaction (5) 
and (7), and the rate constants are 1.04 x 10 cm -molecule -s and 1.79 x 10 s , respectively. 
In Channel B, the rate constants of the abstraction reactions (12) and (14) are 8.70 x 10" 50 and 
9.84 x 10 cm -molecule -s , which are difficult to achieve because the reaction barrier is high. 

For the rate constants over the temperature range of 200-500 K, the Arrhenius equations, 
i.e., k(T) = A exp(—Ea/RT), are shown in Table 3. The pre-exponential factor and the activation energy 
can be obtained from Arrhenius equations. The correlative coefficient R is above 0.997. 

Table 3. The Arrhenius equations for the rate constants k(2)-k(\5) over the temperature 
range of 200-500 K. 



Reaction 


A 


Ea (kJ/mol) 




Arrhenius Equation 


R 2 


(2) 


4.66 x 10 12 


2398.76 


k 


= 4.66 x 10 12 exp(-288.52/7) 


0.997 


(3) 


9.13 x 10~ n 


5564.81 


k- 


= 9.13 x 10" exp(-669.33/7) 


0.999 


(4) 


3.22 x 10 12 


207,450.90 


k 


= 3.22 x 10 12 exp(-24,952/7) 


0.999 


(5) 


1.30 x 10~ 16 


23,352.36 


k = 


= 1.30 x 10~ 16 exp(-2808.80/7) 


0.999 


(6) 


8.05 x 10 10 


118,399.70 


k 


= 8.05 x 10 10 exp(-14,241/7) 


0.999 


(7) 


7.89 x 10 10 


163,511.40 


k 


= 7.89 x 10 10 exp(-19,667/7) 


0.999 
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Table 3. Cont. 



Reaction 


A 




Ea. (kJ/mol) 




Arrhpnins Kniiatinn 


R 1 




6.35 x 10" 


16 


33 661 72 


k ~- 


= 6 35 x If) -16 exr>(-4048 80/71 


0.999 


(9) 

\y) 




12 


1 36 798 60 


k 


= 2 30 x 1 0 12 exnf-1 6 454/71 


0 999 


(10) 


3.18 x 10 


12 


21,194.88 


k 


= 3.18 x 10 12 exp(-2549.30/7) 


0.999 


(11) 


8.98 x 10" 


-15 


11,478.31 


k 


= 8.98 x 10 15 exp(-1380.6/7) 


0.998 


(12) 


7.04 x 10" 


-32 


102,328.70 


k 


= 7.04 x 10 _32 exp(-12,308/7) 


0.998 


(13) 


5.49 x 10" 


-14 


27,305.67 


k 


= 5.49 x 10~ 14 exp(-3284.3/7) 


0.999 


(14) 


1.09 x 10 


16 


205,547.02 


k 


= 1.09 x 10~ 16 exp(-24,723/7) 


0.999 


(15) 


7.12 x 10" 


16 


85,567.69 


k 


= 7.12 x 10 _16 exp(- 10,292/7) 


0.999 



k(T) = A exp(-Ea/RT), where A, pre-exponential factor; Ea, activation energy; R, ideal gas constant 
(R = 8.3 14); T, temperature (K); R 1 , correlation coefficient. 

2.3. QSAR Models 

The quantitative structure-activity relationship is performed to reveal the relationship between the 
structures of PFCAs radicals (~Cs-C2) and the rate constants of decarboxylation reactions at 298.15 K. 
The atom number of PFCAs (2 < n < 7) radicals is drawn in Figure 3. The obtained parameters, such as 
the bond length, molecular mass, volume, dipole, bond energy and net atomic charge, are listed in Table 4. 
Table 5 shows the actual values, predicted values and residual values of the model. The comparison of 
actual values and predicted ones is shown in Figure S2. 

Figure 3. The atom number of PFCAs (2 < n < 7) radicals. 



(2^«^7) 



The independent variables are chosen as follows: the bond length of C1-C2 (i?ci-C2), the molecular 
mass (M), the volume (V), the difference AE between £homo and £lumo, the dipole of molecule, the 
bond energy of C1-C2 (itci-C2), the net atomic charges on atom Ci and C2 (Qci and Qci)- The dependent 
variable is the logarithmic form of rate constants (log k). Then, the genetic function approximation 
(GFA) is adopted to filter the optimum combination of parameters. A five-parameter combination is 
obtained, which includes sufficient information and high reliability in predicting the rate constants of 
decarboxylation reactions. The structure-activity model obtained from the GFA calculation is given 
as follows: 

log k= 361.39 xi? cl _ C2 - 0.03 x F-2.87 x AE + 0.22 x Eci-cz - 27.73 x Q C2 - 534.97 (1) 

The correlative coefficient R 2 and the cross validated R 2 (CV) are above 0.999. It can be seen from 
Table 5 that the relativity between actual and predicted values is excellent. 
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As the structure-activity model shows, i?ci-C2 and £ci-C2 are positively correlated to rate constants, 
while V, AE and Qa are negatively correlated. It is obvious that the rate constant is largely affected by 
i?ci-C2 due to the high factor. When the bond of i?ci-c2 is elongated, the stereo effect is reduced, and the 
rate constant is increased correspondingly. 



Table 4. The parameters used to make quantitative structure-activity relationship analysis. 



PFCAs 
Radicals 


^Cl-C2 

(A) 


M 


V 

(cm 3 /mol) 


AE 
(a.u.) 


Dipole 

(Debye) 


Eqi-ci 
(kcal/mol) 


fici (C) 


fic2 (C) 


C 7 F 15 COO 


1.5433 


412.9659 


166.972 


0.3634 


1.1594 


24.12 


0.659 


0.365 


C 6 F 13 COO 


1.5386 


362.9691 


152.634 


0.3616 


1.1382 


24.22 


0.663 


0.367 


C 5 F„COO 


1.5388 


312.9723 


141.161 


0.2728 


1.2436 


24.15 


0.600 


0.374 


C 4 F 9 COO 


1.5441 


262.9755 


118.218 


0.2533 


1.2094 


24.10 


0.455 


0.457 


C 3 F 7 COO 


1.5394 


212.9787 


83.114 


0.3665 


1.2019 


24.26 


0.472 


0.417 


C2F5COO 


1.5380 


162.9819 


62.998 


0.2863 


1.1352 


23.03 


0.322 


0.418 


CF3COO 


1.5392 


112.9850 


46.759 


0.2894 


1.1828 


19.69 


0.936 


0.419 



Table 5. The actual values, predicted values and residual values of the model (T= 298.15 K). 



PFCAs Radicals 


Lg£ (Actual Values) 


Lg& (Predicted Values) 


Residual Values 


C 7 F 15 COO 


12.2430 


12.2439 


-8.34 x 10" 4 


C 6 F 13 COO 


10.9015 


10.9029 


-1.42 x 10~ 3 


C 5 F„COO 


11.3522 


11.3497 


2.43 x 1(T 3 


C 4 F 9 COO 


11.6284 


11.6286 


-2.02 x 10~ 4 


C 3 F 7 COO 


11.7627 


11.7603 


2.43 x 1(T 3 


C 2 F 5 COO 


11.7340 


11.7373 


-3.28 x 10~ 3 


CF 3 COO 


11.8579 


11.8571 


8.72 x 10~ 4 



3. Experimental Section 

3.1. Geometry Optimization 

The geometrical parameters of reactants, transition states, intermediates, and products are optimized 
at the MPWB1K/6-3 1 + G(d, p) level. The vibrational frequencies have been calculated at the same level 
in order to determine the nature of stationary points. The MPWB1K method is a hybrid density 
functional theory (HDFT) model developed by Truhlar et al. Study shows that MPWB1K gives 
the best results for a combination of thermochemistry, thermochemical kinetics, hydrogen 
bonding and weak interactions, especially for thermochemical kinetics and non-covalent 
interaction [30]. Compared with other conventional methods such as B3LYP and MP2, MPWB1K is 
more accurate and less time-consuming. The 6-31+G(d, p) basis sets are chosen to perform the geometry 
optimization and the 6-311 + G(3df, 2p) basis sets are adopted to calculate the potential energy for this 
medium-scale system after overall consideration of the computational accuracy and time cost [31,32]. 
Each transition state is verified to connect the specific reactants with products by performing an intrinsic 
reaction coordinate (IRC) analysis. All the calculations are performed using the GAUSSIAN 03 
programs [33]. In this study, TS, IM and P represent the transition state, the intermediate and 
product, respectively. 
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3.2. Kinetic Calculation 

Among the minimum energy path, about 40 points near the transition state are selected to perform 
the vibrational frequency calculation, 20 points on the reactant side and 20 points on the product side 
which should represent the shape of the minimum energy path. Based on the information from ab initio 
calculations, including coordinates, gradients, and force constants or Hessian matrix, the rate constants 
with tunneling effects are calculated by the POLYRATE 9.7 program [34]. The canonical variational 
transition state theory (CVT) with small-curvature tunneling (SCT) effect correction is a useful 
method to calculate rate constants [35], which has been successfully applied to lots of research [36-38]. 
The Rice-Ramsperger-Kassel-Marcus (RRKM) theory is used to calculate the rate constants of the 
reactions that have no transition states. Study shows that the rate constants of barrierless reactions 
calculated by RRKM are in good agreement with the experimental observation [39]. 

3.3. QSAR Analysis 

The genetic function approximation (GFA) [40] in the Materials studio package is adopted to 
describe the relationship between the rate constants of decarboxylation reactions and the structures of 
PFCAs radicals (-C8-C2). The parameters obtained from geometry optimization and frequency 
calculation, such as the bond length, molecular mass, volume, atomic net charge and frontier orbital 
energies (Ehomo, ^lumo, A£), are chosen as the independent variables. The dependent variable is the 
logarithmic form of rate constants, log k. The optimal parameter combination can be obtained after analysis. 

4. Conclusions 

In this paper, the rate constants of electrochemical reactions are calculated using RRKM theory and the 
CVT with SCT correction. The structure-activity relationship is analyzed in order to find out the relationship 
between the structures of PFC A radicals (-C2-C8) and the rate constants of decarboxylation reactions. 

(1) The circle of C„F2„+i — > C„-iF 2 „-i is the optimal reaction pathway in the degradation process due 
to the high rate constants. 

(2) The quantitative structure-activity relationship is investigated using the GFA method. 
The structure-activity model has been constructed: i?ci-C2 and E C i-ci are positively correlated to the rate 
constants of decarboxylation reactions. 
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